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Ps| Abstract. This work analyzes the effects of cubic nonlinearities on certain resonant 

r^ scattering anomahes associated with the dissolution of an embedded eigenvalue of a linear 

pvj scattering system. These sharp peak-dip anomalies in the frequency domain are often 

^ called Fano resonances. We study a simple model that incorporates the essential features 

^2 of this kind of resonance. It features a linear scatterer attached to a transmission line 

^ with a point-mass defect and coupled to a nonlinear oscillator. We prove two power laws 

T-j- in the small coupling (7 — t- 0) and small nonlinearity (^ — t- 0) regime. The asymptotic 

relation /x ~ 6*7^ characterizes the emergence of a small frequency interval of triple 

(— I harmonic solutions near the resonant frequency of the oscillator. As the nonlinearity 

Qh grows or the coupling diminishes, this interval widens and, at the relation /u ~ C'7^, 

r"! merges with another evolving frequency interval of triple harmonic solutions that extends 

jrt to infinity. Our model allows rigorous computation of stability in the small // and 7 limit. 

^ In the regime of triple harmonic solutions, those with largest and smallest response of 

I— I the oscillator are linearly stable and the solution with intermediate response is unstable. 
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rt- nance, bist ability, 
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^ 1 Introduction of the nonlinear model 

^ The interaction between a resonant scatterer and an extended system that admits a spectral con- 

tinuum of states is a fundamental problem in classical and quantum systems. A variety of simple 
models have been devised to elucidate this interaction, and they have the advantage of providing 

C^ clean mathematical treatments of specific phenomena. H. Lamb [8j observed in 1900 that, if a 

simple harmonic oscillator is attached to a point on a string, the loss of energy by radiation into 
the string effectively results in the usual damped oscillator. Komech [7] extended the analysis to 
a general simple nonlinear oscillator and proved that finite-energy solutions tend to an equilibrium 
state of the oscillator and in fact that transitions between any two equilibrium states are possi- 
ble. In the frequency domain, resonance effects of an oscillator on extended time-harmonic states 
were treated by Fano [1] (1961), in order to explain peak-dip anomalies (the "Fano line shape") 
observed in the scattering of electrons by the noble gases. These sharp resonances are a result of 
a weak coupling of a bound state to a continuum of extended states. The frequency of the bound 
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state is realized as an eigenvalue of the equations of the extended system, embedded in the con- 
tinuous spectrum, and this eigenvalue is unstable with respect to perturbations of the system. An 
analogous phenomenon occurs in the scattering of an EM plane wave by a photonic crystal slab. 
Sharp peak-dip anomalies in the transmission of energy across the scatterer (similar to those in 
Fig. p] and often referred to as resonances of Fano type) are due to the resonant excitation, by the 
incident wave, of the state that is localized in the slab [H [2l [16] . These resonances can be analyzed 
rigorously by means of a complex-analytic perturbation theory of the scattering problem about the 
parameters of the bound state, and one obtains asymptotic formulae that reveal fine details of the 
anomalies |14l I17j . The analysis is applicable quite generally to scattering problems that admit 
unstable bound states, including continuous as well as lattice models [TSl fTS]. 

Kerr (cubic) nonlinearity in models of resonant harmonic scattering have been investigated by 
several authors, mostly numerically. Its effects on resonance in dielectric slabs is important for 
applications exploiting tunable bistability [9] . The Fano- Anderson model of a resonator coupled to 
a single chain of "atoms" with nearest-neighbor interactions exhibits a sharp dip in the transmission 
coefficient in the middle of the spectrum. When Kerr nonlinearity is introduced into the resonator, 
a stable scattering state bifurcates as the nonlinearity reaches a critical value, producing multiple 
scattering solutions and bistability in an interval of resonant frequencies [12]. Discrete nonlinear 
chains have been proposed for modeling photonic crystal waveguides [1U\ [TT] . When, instead of a 
single atomic chain, two atomic chains are coupled together, one can construct embedded trapped 
modes exhibiting the Fano- type anomaly observed in photonic systems, and numerical computations 
show that Kerr nonlinearity causes intricate multi- valued transmission coefficients |15j . 

In order to understand the fundamental effects that nonlinearity produces upon resonant scat- 
tering through exact explicit formulae, we introduce a dynamical system consisting of two systems 
coupled together (Fig. [T]): (1) a transmission line of Schrodinger type with a point-mass defect 
acting as a non-resonant scatterer and (2) a nonlinear resonator coupled to the defect on the line. 
When the nonlinearity is set to zero, the model incorporates the essential features of these linear 
photonic systems described above inasmuch as the resonant phenomena are concerned. The model 
elucidates fundamental nonlinear effects that coincide with asymptotic relations between the pa- 
rameters of coupling and nonlinearity. A discussion and comparison with the models of Lamb and 
Komech as well as the Duffing oscillator is offered in the final section [5] The reader may decide to 
browse that discussion before embarking on the analysis of our model in sections [2j|4] 

The transmission line with a point mass at x = is the system 

iut{x,t) = -huxxix,t), X / 0, 
mi(0,t) = -r(u^(0+,t) -n^.(0-,t)), 

in which /i > carries units of area per time and r > carries units of length per time. Selecting 
T = h/r"^ and X = hjr as units of time and space nondimensionalizes the equation and removes 
h and r. Thus the nondimensional frequency w = 1 corresponds to the dimensional frequency 
27:7"^ /h. The resonator is a nonlinear harmonic oscillator with cubic, or Kerr, nonlinearity, that 
obeys the equation in nondimensional form 

iz{t) = Eoz{t) + X\z(t)\'^z{t). 

The nondimensional characteristic frequency Eq is relative to the fixed frequency unit 27r/T, and 
A is a nondimensionalization of a parameter carrying units of frequency per square units of z (the 
values u and z of the fields are assumed to be nondimensionalized relative to a common unit). 
Neither £"0 nor A can be eliminated from the equation. 



These two systems are then coupled by a parameter 7, 

iut{x,t) = —Uxx{x,t), X ^ 0, (transmission line) 

iy{t) = 'yz{t) — {ux{0~^,t) — Ux{0~ ,t)), y{t) = u{0,t) (point defect on the line) (1-1) 

iz{t) = Eoz{t) +jy{t) + X\z{t)\'^z{t). (resonator) 

Nonlinearity is measured by the composite parameter 

where J is the amplitude of a monochromatic field emanating from a source at —00 (see Fig. [l] or 



(2.2)). In its dimensional form, /i and 7 are frequencies, with /x being amplitude-dependent and 
depending quadratically on the strength of the incident field. 

The resonant amplification of the oscillator in the linear system (A = 0), accompanied by sharp 
transmission anomalies near the resonant frequency Eq, as described above, is portrayed in Fig. [2j 
As the coupling 7 of the oscillator to the transmission line vanishes, the resonant amplification 
becomes unbounded. In the limit 7 — )• 0, the motion of the oscillator becomes a bound state for 
the decoupled system, whose frequency is embedded in the continuous spectrum arising from the 
transmission line. 

The nonlinear system (1.1) admits time-periodic solutions, which the Kerr form A|zpz of the 



nonlinearity allows to be monochromatic. Even a small nonlinearity has a pronounced effect on 
the resonance, due to the high-amplitude fields produced. Fig. |3] shows how multiple scattering 
solutions (nonuniqueness of the scattering problem) appear near the resonant frequency and spread 
to higher frequencies as the nonlinearity parameter is increased. Similarly, a new branch of solutions 
emerges from the infinite frequency limit and spreads to lower frequencies. 

The present work proves that the bifurcations occur as portrayed in Fig. |4] as the nonlinearity 
increases from the value zero. The analysis reveals characteristic power laws between the coupling 
parameter 7 and the composite parameter of nonlinearity ^ in the asymptotic regime of 7 — t- and 
/i — 7- (see Table 111 and Fig. Q. The critical relation 7^/// ~ Const, marks the onset of a narrow 
frequency interval [cji,a;2] of triple solutions near the resonant frequency. For each frequency above 
this interval and below a high frequency cja ~ C'y~'^, the system admits a unique harmonic solution. 
As fi increases relative to 7, the intervals [u;i,W2] and [a;3,c«) of triple solutions widen until they 
merge, that is a;2 = W3 ~ C, at the critical power law j'^/fi ~ Const. For any fixed, small /i 
and 7, there are at most three frequencies that separate intervals of unicity from intervals of triple 
solutions; we call these transition frequencies. Understanding of the dependence of the transition 
frequencies cji, LJ2, and ws on the parameters 7 and /j, can be illustrated graphically through the 
intersection points of two curves in the plane whose coordinates are reparameterizations of the 
frequency of a harmonic excitation and the "response" of the resonator (§3|. 

The stability of multiple scattering solutions is important in applications involving bistable 
optical transmission [31 IHl [T^ I19| . It is commonly understood that, at frequencies for which a 
nonlinear resonant system admits a unique periodic solution (for a given system and source-field 
amplitude J), this field is stable under perturbations. At frequencies for which there are exactly 
three periodic solutions, there should be bistability, that is, two solutions are stable and one unsta- 
ble, the unstable one being that with intermediate amplitude (see, e.g., the references just cited). 
The analysis in this section yields precise statements about this kind of bistable behavior in our 
model of nonlinear resonance (Theorem [2]) . 
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Figure 1: The model for nonlinear resonant scattering. The transmission line models the ambient space, 
the mass attached to the line models the (non-resonant) scatterer, and the nonlinear oscillator attached to 
the mass models the bound state, which is responsible for resonance. The Kerr nonlinearity is confined to 



the oscillator. The elements of a harmonic scattering field at frequency uj correspond to the form (2.2 1 




Figure 2: Transmission versus frequency for the linear system (A = 0) and normalized incident amplitude 
J = 1 for three values of the coupling parameter 7: From left to right, 7 = 0.0,0.2,0.7. The characteristic 
frequency of the oscillator is cu — Eq ^ A. 
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Table 1: Asymptotics of the bifurcations of the transition frequencies for small 7 and /_*. The bifurcation 
diagram with /i as the bifurcation parameter and fixed 7 is shown in Fig. |4] All constants depend only 
on Eq and are computed explicitly in the text. The three relations in the first row are given in equations 
(3.24|3.25|3.26); the relations in the second row are given in equations (3.31|3.34|3.30|. 



2 Harmonic solutions of the nonlinear model 



System (1.1) admits harmonic scattering solutions, in which a field Je*' ^^' oscillating at fre- 



'ikx\ —iict 



quency uj in the string is incident on the resonator from the left (Fig. [T]), 

u{x,t) = (Je*'=^ + ^e- 
u{x,t) = Be'^''e-^'^\ 
y(t) = Be-''^\ 
z{t) 



X <0, 
X > 0, 



(2.2) 
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where y{t) = u(0, t) represents the state of the defect in the transmission line. The wave number 
/c > in the string and the frequency uj are subject to the dispersion relation for the free string 
UJ = k'^. The continuity of u provides the relation 

J + A = B, 



and (1.1) reduces to the algebraic system 

{2ik-k^)B + -fZ = 2ikJ, 

-fB + {Eo - k^)Z + X\Z\^Z = , 

which can be written as a cubic equation for Z/J and linear a relation between B/J and Z/J : 
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(2.3) 



Evidently the dependence on A and J is only through the composite parameter 

^l := XJ^. 

In Figures. [3]and[5j one observes two intervals of the w-line in which the harmonic scattering 
problem has three solutions and two intervals in which it has one solution. These intervals are 



separated by three transition frequencies kf = uJi{'y, fj,,EQ 



1,2,3: 



(0, UJl) U (a;2, W3) one solution, 
(ci;i, ^2) U (^3,00) three solutions. 



(2.4) 
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Figure 3: The transmission coefficient T := \B/ J\'^ of tfie harmonic scattering solutions ( |2.2[ ) of the string- 
oscillator system ( |1.1[ ) us. the wavenumber fc = -^07 . The graphs show progression of values of the composite 
nonlinearity parameter ^ = AJ^ (from top to bottom, /^t = 0, 0.001, 0.0035, 0.003856, 0.00395), all other 
parameters being fixed [Eq = 4, 7 = 0.3). The top graph shows the linear case, in which each frequency 
admits a single harmonic solution. The sharp anomaly occurs near the characteristic frequency i?o = 4 
of the oscillator. For nonzero jjl, there is a half-infinite frequency interval {lo^^oo) of triple solutions. As 
/.t increases (or as 7 decreases), a narrow frequency interval {ijji,ui2) of triple solutions emerges; this is the 
U)i-uj2 bifurcation. Meanwhile, lo^ decreases, eventually merging with 1^2 and eliminating the interval {u2, W3) 
of unique solutions; this is the a;2-W3 bifurcation. A diagram of the bifurcations for this progression is shown 
in Fig. H (bottom). 
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Figure 4: Top: Transition frequencies wi, a;2, and ws, vs. /i, at which the number of harmonic scattering 



solutions of the form (2.2 1 changes from 1 to 3 or vice- versa, according to whether the cubic polynomial (2.8) 



for the response has one or three real roots. The unshaded region indicates those {uj,h) pairs that admit 
only one real root, whereas pairs in the lightly shaded region admit three real roots. Bifurcations occur at 
the points B12 (the cuspidal U!i-u}2 bifurcation in the text) and B23 {uj2-^3 bifurcation) as fj, varies. The 
coordinates of the two bifurcations are labeled with their asymptotic values for small 7 and fj,, according to 
Table [1] Bottom: This is a numerical computation for Eq = 4 and 7 = 0.3 that corresponds to Fig.|3] 



When /i is decreased or 7 increased, the points wi and L02 approach each other and are annihilated, 
and there remains a single interval of one solution and one of three solutions, separated by the 
point UJ3. We will call this threshold the 001-002 bifurcation. On the other hand, increasing p. or 
decreasing 7 brings the points 002 and 00^ together, resulting in a single interval [ci;i, 00) of multiple 
solutions. We call the threshold that occurs when they are annihilated the 002-00^ bifurcation. 

Each of these bifurcations is analyzed in detail in the following section, and Table[T]is a summary 
of the results. Asymptotically, as /i and 7 vanish, the bifurcations are characterized by specific power 
laws. The a;i-a;2 bifurcation occurs at relatively weak nonlinearity, p ~ C't^j when its effects on the 
resonance near Eq first appear. The 002-00^ bifurcation occurs at much higher values of nonlinearity 
relative to the coupling parameter. These constants as well as those in the asymptotics for the 
frequencies ooi are computed explicitly in terms of the resonant frequency Eq alone. 

Analysis of the number of scattering solutions as a function of frequency is facilitated by a 
polynomial equation for the square of the amplitude enhancement TZ = \Z/ J\'^ in the resonator. 



which we call the response. The first equation of (2.3) yields 



/ ^2 \2 4 2 

7^U7^+^— -(w-^o) + , ,,0 ^ TT = 0- (2-5) 

V a; + 4 ^ 7 oo{u) + AY a; + 4 

The roots TZ of this equation, which are necessarily positive, are in one-to-one correspondence with 



the solutions {B,Z) of (2.3). This equation can have multiple roots only if 

{00 — Eq){oo -|- 4) — 7 > 0, (necessary for multiple solutions) (2-6) 

in particular, if a; > Eq^ that is, if 00 exceeds the frequency of the free resonator. The change of 
variable 

n=-(uj-Eo- ^— ) (2.7) 



transforms equation (2.5) into 

/(P):=P(P-l)2+aP-/3 = 0, (2. 



in which 

1 / 272 



a = a(a;;7,£'o) := 



/3 = P{u;;-f,p,Eo) : 



00 \{uo - Eq){oo + A) - -f"^ 
Aj'^nioo + 4)2 



(2.9) 



The quantities 7^, a, and /3 can be written more compactly using the notation 



p := oj + 4 


r^ 00 


as w — )■ 00, 


a := 00 — Eq 


— 72/p ~ 00 


as w — )> 00. 


a 

n = -p, 


4^4 

a= , 
ooa^p'^ 


a-^p 



If one fixes the parameters p, 7, and Eq of the dynamical system, one can consider the roots of 
the polynomial f{P) as functions of frequency co. These roots correspond to all possible responses 7^ 
associated with harmonic scattering solutions at 00. Transition frequencies separating unicity of the 
response from multiple responses occur when f{P) has a double root. Thus analysis of transition 
frequencies is tantamount to analysis of the parameters for which f{P) has a double root. 
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Figure 5: These graphs show the harmonic scattering states (2.2 1 of the string-oscillator system (1.1 1 vs. 
the wavenumber k = ^JZj. Plotted from top to bottom are (i) the logarithm log|Z| of the amplitude of the 
oscillator's motion z{t) — Ze~'"*; (ii) a rescaling P of the response 7?. = |^/Jp defined by (2.7) (J is the 



amplitude of the incident field) with detail near the resonant frequency depicted to the right; (iii) the phase 
arg(Z) of the oscillator; and (iv) the transmission coefficient T— jB/ Jp. For the parameter values i?o =4.0, 
7 = 0.3, and /i = 0.0035 chosen here, there are two frequency intervals (wi,a;2) and (^3,00), delineated by the 
vertical dotted lines, in which the system has three solutions; at frequencies in the complementary intervals 
(0,aji) and (aj2,a;3), the system has only one solution. The two high-amplitude oscillations have almost 
identical amplitudes, as shown in the top graph, but they differ by a phase that tends to tt as a; — )■ cxi, as 
seen in the third graph. A graph of the transition frequencies vs. [i is shown in Fig. |4] 



3 Multiple harmonic solutions and their bifurcations 

We have seen that certain transition frequencies uji(iJ,,^;EQ) separate intervals of unique response 
to harmonic forcing by the source field Je"*'^* from intervals of triple response. The diagram in 
Fig. |4] shows how the transition frequencies merge or separate as a function of the nonlinearity in 
the asymptotic regime of 7 — t- and /i — t- 0. This section is devoted to proving the asymptotic 
power laws in that diagram. The main results are stated in the following theorem. 

Theorem 3.1. The transition frequencies cji that separate frequency intervals of unique scattering 
solutions from intervals of multiple solutions vary with the nonlinearity parameter // according to 
the diagram in Fig. [^ if the coupling parameter 7 is sufficiently small. More precisely, 

1. There are at most three transition frequencies if 'j is sufficiently small. 

2. There is a point B12 = (wi2,/^*) and positive numbers C12 and C'12 with 

— ~ C12 and uj\2 — £'0 ~ C'[97^ (7 — ^ 0) 
/"* 

such that, if [I < //=„, there is a single transition frequency uj^ and if fi > jjl^ hut /j, is not too large, 
there are three transition frequencies with oji < U2 < UJ3. ^s /i — )• /i=K from above, ui — )• UJ12 and 
L02 — )■ ooi2- Moreover when ^ = ^^, there are two transition frequencies UJ12 = loi = i02 < ^"i and 
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Cv2 

1' 



(7-^0). {UJI=UJ2) 



The numbers C12 and C'12 depend only on Eq and are given by (3.24) and (3.25). 
3. There is a point B23 



{u}23, n*) and positive numbers C23 and C 23 with 

„2 



1_ 



C23 and 6^23 ~ C23 (7 -^ 0) 



such that, if fi > fi* , there is a single transition frequency wi and if fi < fi* but /x is not too small, 
there are three transition frequencies with wi < L1J2 < W3. As fj, ^ fi* from below, uj2 — ^ W23 and 
UJ3 — )• (jj23 . Moreover, when /x = /x* , there are two transition frequencies uji < 0J2 = ijJ3 = a;23 and 



wi — £"0 



C, 



1/3 
23 



7^/3 (7^0). 



The numbers C23 and C23 depend only on Eq and are given by (3.31) and (3.30). 



The transition frequencies are characterized by the property that there exists a real number P 



such that both (2.8) and its derivative with respect to P vanish. The system of the two conditions 



is algebraically equivalent to the pair 



a(a;;7,£o) 
/3(w;7,/x,£o) 



-3(P-|)(P-1), 

-2P2(p_l). 



(conditions for transition frequencies) (3.10) 



It is convenient to work with the roots of ap = {u — Eq){uj + ^) — 7^ = {uj — Ei){(jJ + c), which 



10 



are small perturbations of Eq and —4 as 7 — )• 0, 



(uj - Ei){uj + c) = {uj - Eo)iuj + 4) - -f"^ , 

El = Eq + e, 

c = 4 + e, 



7 



+ 0(7' 



a = 4 + El 



Eo + 4 

4 + Eo + €. 



Because of inequality (2.6), multiple solutions are possible only for uj > Ei, and 
introduce the variables 

V = bJ — El, 

q = 3P-2. 



(3.11) 
(3.12) 
(3.13) 

(3.14) 

(3.15) 

we therefore 

(3.16) 
(3.17) 
(3.18) 



The algebraic system (3.10) in P and oj can be rewritten as a system in q and r: 

^5 



(l-^)(l + g) = 3a = 127^ 



T 



[l + EiT){l + {a + e)Ty 



2 _ 27 547V-r^(l + aT)2 



{l-q){2 + qr=-^(3- ,,^, ^ , ,s ■ 
2 (l + (a + e)r)'^ 

Dividing the second by the first and retaining the first equation yields the equivalent 



pair 



f {l-q){l + q) -- 

1 i^ + Q? 
[ l + q 


_ io.,4 r^ _ 127^ 




(a) 


^' {1 + Eit){1 + {a + e)Ty u^u + Ei){u + {a + e))^ ' 

9 fi {1 + ar)2(l + Eit) _ 9 fi {u + af{u + Ei) 
272 r(l + (a + e)r) 272 i/(i/ + (a + e)) 



(3.19) 



3.1 Graphical depiction of transition frequencies and their bifurcations 

These two algebraic relations between q and r are shown in Figs. [7] and |8j the symmetric one is 
(3.19i). They provide a transparent graphical means of analyzing the transition frequencies and 



their bifurcations. The r-values of the points of intersection between the two relations determine 
these frequencies through u: = Ei+t~^ . It is visually clear that there are at most three intersection 
points, and we give a proof of this in section |3.5[ Fig. [7] shows the evolution of the relations as /i 
increases, for a fixed value of 7. Initially, there is a single intersection, corresponding to ujz- At a 
critical value of /i, another intersection appears at the top, which then splits into two intersections 
corresponding to toi {q < 0) and UJ2 {q > 0). This is the U1-UJ2 bifurcation. After a mountain-pass 
mutation of relation (3T9p), the intersection points corresponding to i02 and uj^ approach each 
other, fuse together, and then disappear on the lower right half of the relation ( 3.19^ ) in the uj2-(^3 
bifurcation. 

When 7 tends to zero, the peak of relation ( 3.19^ ) grows without bound. But the frequency of 
the u;2-'^3 bifurcation remains of order 1 and (3.19 3,) appears as two practically vertical lines, one 
at q = —1 and one at q = 1. This regime is depicted in Fig. [8 

The graphical realizations of relations (3.19) in Fig.sITland qare obtained as follows. In ( 3.19^ ), 
the function of r maps the positive real line onto itself in a strictly increasing manner. Thus r 
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Figure 6: The concave functions fi{q) and /2(t) in the relation (3.19 d). The function /i(g) does not depend 



on the parameters of the system, whereas /2(t) does. When the minimal value of /2 is less than that of /i 
(equal to 4), the relation fi{q) = /2(t) possesses two components that are functions of q, with minimal and 
maximal r- values equal to t^ and r+, as in the first three graphs of Fig. [tI When min/2(r) > 4, the two 
components are functions of r, with minimal and maximal r-values equal to g_ and g-|_ as in the last three 
graphs of Fig. [7j 




Figure 7: The relations (3.19) in the ( g, r) plane, whose intersections give the transition frequencies through 
uj — El + r~^. The symmetric one is (3.19 1). The values of 7 and Eq are fixed, and the graphs from left to 
right show the evolution of the relation (3T9b) as fi increases. The two bifurcations occur at the tangential 
intersections, as in the second {uji-ll)2 bifurcation) and fifth (UJ2-UJ3 bifurcation) graphs. 



12 



is implicitly a function of g G (0, 1) with infinite slope at q = ±1, and we denote this symmetric 
function by r = T{q). 

Relation (3.19[3) can be understood by placing the graphs of the functions 



/2(r) 



(2 + g)^ 

9 /i (l + ar)^(l + giT) 
272 r(l + (a + e)r) ' 



(3.20) 
(3.21) 



side by side, as shown in Figure |6j The function /i is convex and tends to infinity as <; — )• — 1 or 
q — )• 00; its minimal value of 4 is achieved at g = 0. Likewise, /i is convex and tends to infinity 
as r — )• or T — )• 00, and thus it has a minimal value, say ttiq. If niQ < 4, then ( 3.19) 3), or 
fi{q) = /2('7")) has two components, each of which is the graph of a function of q, one with unique 
local maximum r_ and the other with unique local minimum r+ , with r_ < T4. , both of which are 
achieved at q = 0. If mo > 4, then each of two components is the graph of a function of r, the 
lower having maximal value q- and the upper having minimal value g+, with g_ < < g+. At the 
transition from one regime to the other, when mo = 4, the relation consists of two curves crossing 
tangentially at g = 0. 

If T is viewed as a multi- valued function of q, then the upper branch of ( 3.19[ 3) is decreasing 
for g < and increasing for q > 0. Since (3.19i) has the opposite behavior, it intersects the upper 
branch of ( 3.19| 3) either not at all, exactly once at q = 0, or at two points, one with q < and 
one with q > 0. One can see that the r- value of the former is larger by the observation that 
fi{~l) ~ fiio) = 2g^/(l — g2) > 0, for q G (0, 1). It is clear from the graphs that the lower branch 



intersects (3.19i) at most once, and only for g > 0; a rigorous statement is proved in section 3.5 



3.2 The UJ1-UJ2 bifurcation 

We first analyze the bifurcation occurring at the point B12 of Fig. [4j at which a narrow frequency 
interval of multi-valued harmonic solutions is born as fi increases across a threshold value (with 
7 fixed). We prove the power law 7^ ~ Ci2^ and the asymptotics of the transition frequencies 
at this bifurcation. The sharp feature in the bifurcation at the point B12 is in fact cuspidal with 
a;2 — wi ~ C{^ — jioY'"^, meaning that the {(jJ\,0J2) interval opens slowly. 

The maximal value of the function r = T{q) is attained at q' = 0, and we denote it by 

T* = v~^ = max T{q) = T(0). 
96(0,1) 



Setting g = in (3.19) gives 

As 7 



127^ = ul{u^ + Ei){v^ + a + ef. 
with v^ > 0, we have v^ — )■ 0, z^* + £'1 — >■ £'0, and v + a + e 

1274 



Eq + 4, and thus 



from which we obtain 



n 



-^Eo{Eo + ^f (7^0) 
^o{Eo + 4) _ C, 



2^/372 



r 



(7 ^ 0). 



(3.22) 



The pair of frequencies oji and uj2 is born (or annihilated) when the upper branch of the second 



relation in (3.19) intersects the first relation r = T{q) at its peak, that is, when r+ = r,,,, as shown 



in the middle top graph of Figure ^ 
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Figure 8: The relations ( 3.19^ b) for small 7 and /x. The intersections of the two relations correspond the 
the transition frequencies. Here, 7 = 0.1 and Eq = 4, and fi increases through values ranging from the loi-(jJ2 
bifurcation at ^i -- 7^/Ci2 « 74/I66 w 6.0e-7, through the crossing at ^i -- I'^/C w 72/26.2 w 0.00038122, 
to the wi-wa bifurcation at ^ - 7VC23 ~ 7^/23.3 « 0.000429. 



The numbers t± are the solutions of the second relation of (3.19) with g = 0. This equation 
can be written as 



r?-— r = (l + aT)(l + SiT) 

9 jJL 



in which 



1 + or + er r 

V ■= — ^- = 1 + 6: 



1 + ar ' 1 -|- ar ' 

so that 7/ — )• 1 as 7 — )• 0. Thus we obtain 



T 1 

0< < 



1+ar Eo + 4' 



aEir'^ + (a + Ei-7]-^]T + l = 0. 
9 /U 



(3.23) 



Setting r = r^K in this equation, with 7 — )• 0, the asymptotic relation r* ~ C^:/^'^ yields the balance 
of two terms, 



'Vn — TT ~ "-^1t:J (wi = a;2,7 -^ 0), 



9 /U 7^ 7' 

which results in the asymptotic power law for the uji-uj2 bifurcation 



fj, lb 



E'^^\Eo + if=:Ci2 ioJi=U2,7^0). 



(3.24) 



Let us denote by UJ12 = uji = UJ2 the frequency at which this bifurcation takes place. Equation 
(3.22) and the definitions of z^ and r give uji2—EQ—e ~ 7^/C=k, and this together with e ~ 7^/(4-|-i?o) 
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yields 



"^12 — Eq 



7 



C 



+ 



12 



1 

4 + K 



7 



2^3 + ve;^ 

\/^(£^o + 4) 



7'Gl2 (7^0). 



(3.25) 



As this bifurcation takes place, the third transition frequency 003 is very large, and one can 
compute its asymptotic value by finding the intersection between r = T{q) and the lower branch 



of the second relation of (3.19). Let us denote this intersection by (^3, T3). The maximal r-value of 
the lower branch is t_, which satisfies (3.23), and is seen to be of order 0(7^), and thus the first 



equation of (3.19) gives q = 1 + 0(7 ). Inserting this into the second of (3.19), gives 

v2 



VT3 



T 
C 



12 



■(l + ar3)(l + SiT3) 



in which r; — )• 1 as 7 — )• 0, or 



1 + 



E^ 



^12 \ 2 

-^^■q T + aEiT 

7 / 

-1 



0, 



which has a solution t^ ~ 7 /Ci2. Finally, u^ = t^ + Ei, which yields 



i03 



C12 



{UJI =U2, 7 ^ 0). 



(3.26) 



The response TZi2 of the field at the frequency CJ12 tends to infinity as 7 ^. This can be seen 
by inserting the asymptotic expressions ( 3.24|3.25 ) into (2.7) with q = 0: 

7^l2 = ( wi2 -Eo ^ ) ^ ~ ^7^0(^0 + 4). 

V UJ12 + 4y 3/x 7^4 

Similarly, the response of the field that is created or annihilated at the transition frequency oj^, 



that is, the field corresponding to the double root of (|2.8|), is found to be 

2Cl2 



7^3 



2 C\2 

3 /u2 "" 3 78 



{UJI =V02,1^ 0). 



Let us why the UJ1-UJ2 bifurcation is a cusp. It occurs when a convex function and a concave 
function intersect tangentially at their extreme values, as seen the second graph of the sequences in 
FiguresMand^ The concave function is symmetric and the convex one is not. Up to order 0{q^), 
these functions can be represented by 

-.2 



T = To-cq- , r > 0, , (3.27) 

r = To-r6 + aq^ + bq^, a > , (3.28) 

in which e is a rescaling of /i. After making the substitutions r6/{a + c) 1— t- 5 and d = h/{a + c), the 
intersection points (g, r) of these two relations satisfy 

5 = q^{l + dq). 

The small solutions of this equation have an expansion in powers of yo , 

gi,2 = ±51/2 -^5 ±1^253/2 + ...^ 

and the corresponding r-coordinates are 

The difference of these is 

Tl -T2 = 2d6^^^. 

Seeing that r = l/{uj — Ei) and 6 \s a, rescaling of /i, the difference W2 — ^i is of order (// — ^qY'"^. 
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3.3 Between bifurcations 

The structural morphosis of the relation fi{q) = f2{'T) occurs when the minima of /i and /2 are 
equal and is characterized by the crossing of two curves at (0, tq) as depicted in Figs. |6] and [SJ The 
number tq is where /2 attains its minimum value, say /2(to) = rriQ. By writing 



ft ^ (l + (J^o + 4)r)(l + J^oT) 



1 + ar 



l + (a + e)r' 



in which ^^(r) — 1 and ??'(t) are both 0(7^) uniformly in r, one finds that 



1 



To 



^Ea[Eo + 4) 



(7^0) 



and that 



mo 



:= /2(ro) ~ 94 (2 + ^0 + \/^o(^o + 4)) (7^0). 



The crossing occurs when g = and r = tq simultaneously in the relation fi{q) = /2(''"), that is, 
when 4 = ttt-q, which yields the asymptotic power law 






(2 + ^0 + VEo{Eo + 4)) (at crossing, 7^0). 



When mo < 4, we have t_ < tq < t+, where r-j- are the extremes of the two branches of the 
relation fi{q) = /2(t). Asymptotically, 



T_ < 



1 



^Eo{Eo + 4) 
When mo > 4, the extremes q± are defined through 

i + q± 



+ 0(7) < r+ (7^0). 



mo. 



Since — 1 < g_ < 0, we have 1 < (2 + g_)^ < 4 and therefore 



1 + 



< 



7 



/i \9(2 + ^o + \/^o(^o + 4)) 



+ 0(7) < 4(1 + g_). 



The point (g_,ro) is the rightmost point on the left branch of the relation fi{q) = /2(t), and 
from the diagrams in Fig. ^ this point is evidently to the right of the graph of r = T{q), as the 
latter is practically a vertical line at q = —1 when r is of order 1. The value of q that satisfies 



T{q) = To is obtained by setting r = ro in the first equation of (3.19): 

. / 12 



(1-'Z)(1 



7 



p3/2 



E'o'^Eo + 4)2(/Eb + vw+^y 



+ 0(7) 



Thus q = —1 + 0(7^) and we obtain 



1 + g < 7 



p3/2 



E'^^\Eo + 4y{^/Eo + VEo + 4.y 



+ 0(7) 



(Tiq) = TO, 7 ^ 0). 



(3.29) 
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3.4 The UJ2-CO3 bifurcation 



Because of (3.29) and the symmetry of T, the two intersections for q > merge when g+ = 1+0(7^) 
and T ^ tq. This gives the frequency of the uj2-u}3 bifurcation 



W23 



Eo + VEoiEo + i) =: a, 



23 



(7^0), 



and, setting (2 + q-|-)^/(l + q+) = mo, the asymptotic power law relating 7 to /i, 



^ 



2 (2 + ^0 + VMEo + 4)) =: G 



•'23 



(a;2 = a;3,7-;> 0). 



(3.30) 



(3.31) 



The response at this bifurcation is obtained from (|2.7|), 



7^ 



23 



^Eo{Eo + 4) _ 2(2 + Eo)/Bop(^T4) + 2^0(^0 + 4) 



/^ 



7^ 



As this bifurcation takes place, the intersection that determines wi has g — )• — 1 as 7 — >• 0. This 
is because, with j'^/fj. ~ C23, relation (3.19 a) becomes stationary as 7 — )■ 0, whereas relation (3.19i) 
contains the scaling factor 7^ on the right-hand side. Using 7 — )• and g — )• — 1, together with 



7 //.i ~ C23, system (3.19) gives 



1 



i + q 



9 (l + ^or)(l + (4 + ^o)r) 
2C23 

67^ 



r 

^5 



which together yield 



Thus r — )• 00, and we obtain 



27 



-7 



(l + ^or)(l + (4 + ^o)T)2^ 



1. 



(3.32) 
(3.33) 



C23 ' (1 + (4 + Eq)t) 



3 ' 



Using now uji — Eq + C(7^) = t ^, we obtain 



oji — £^0 



G 



1/3 
23 



v4/3 



(w2 =a;3, 7 -^ 0). 



(3.34) 



The response at the transition frequency wi, corresponding to the double root of (2.8), is TZi 
{loi - Eo + 0{'j'^))/{3ij.), which results in 



7^1 



2 a 



2/3 

23 



3 72/3 



3.5 Proof of at most three transition frequencies 



This section contains the proof of part (1) of Theorem 3.1 It can be reformulated as follows. 



Proposition 1. The system (3.19), or, equivalently, the system (3.10), has at most two solutions 



with < q < 1 (t^ < P < I), and if jjl and 7 are sufficiently small, then the system has at most one 
solution with — l<g<0 (^ < P < ^). 
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Using the relations 3a 
equivalently as 



I- q^ and 27,3/2 = A - 2,q^ - q^ , the system (|3.19|) can be written 



2q 



127^ 






(3.35) 



1 + 



9 



T i^iiy + ia + e)) 

and the first shows that \q\ < 1. 

First we deal with < q < 1. Setting G and H equal to the right-hand-sides of the first and 



second equations of (3.35) gives G as an increasing function of v that maps (0, 00) onto [— c«,l) 
and thus if is a well-defined function of G G (—00, 1). One computes that 



d^H 


UG\-' 


'd^H 


d^GdH 


(dG\-^' 


dG^ 


\du) 


dl/2 


du'^ du 


Kdi^J 



and that d^H/dv"^ > 0, d^G/dv"^ < 0, and dG/du > (for i^ > 0). Thus, whenever dH/du > 0, 
d'^H/dG^ > also. Since 

dH__dH fdCy 

dG ~ 



-1 



du \dv J 



, dG 
and — — > 0, 
dv 



we see that d'^H/dG'^ > whenever dH/dG > 0. 

Let vq be the positive value of u that corresponds to G = 0, so that [z^OjCo) maps onto the 
G-interval [0, 1). Since H is a convex function of v that tends to infinity as 1/ — )• or z^ — )■ 00, it 
has a unique local minimum on [fOjCo) (possibly at vq), and thus H has a unique local minimum 
as a function of G (possibly at G = 0); denote this function by if = 7i{G). 

The expressions G = q^ and H = 2q'^/{l + (7) for g e [0, 1) constitute a parameterization of the 
relation H = 2G/(1 + vG) for G G [0, 1), in which H is an increasing concave function of G. It 
remains to count the number of intersections between H = ?^(G) and H = 2G/(1 + vG) on the 
G-interval [0, 1). Since 'H{G) has a unique local minimum and is convex whenever it is increasing, 
it intersects H = 2G/(1 + vG) no more than twice. 

Now let us consider — 1 < g < 0. Let r = F{q) be defined through the relation fi{q) = f2{T) for 
q G (0,g_) (g_ < 0) and r < tq. We shall show that, for 7 sufficiently small and fi = const. 7, we 
have F"{q) > 0, 1 + g_ > const. 7, and T{q) = tq =^ 1 + q = 0(7^). From this, it follows that the 
relation fi{q) = f2{'T) intersects r = T{q) exactly once for g < and these values of 7 and ^. The 
result is extended to /x < 0(7) by the observation that, as /x decreases, q- increases (until it reaches 
0) and F{q) decreases, which disallows the occurrence of any additional intersections (Fig. |9]). 

We have already found that tq ~ l/ySo(^o + 4), that 1 + q- > const. 7^/;U, and that tq = 



Up) 



1 + q- < const. 7 as 7 — )■ 0. Define 



51(g) 



1 



/i('Z)' 



52 (r) 



{1 + Eot){1 + {Eq + A)t) 



so that F{q) satisfies the asymptotic relation 



9i{q)---92{F{q)) (7^0). 

1^ 



The second derivative of F satisfies asymptotically 

1 



F"{q) 



9'2ir) 



III N /^ 9i{q) II, X 



ir = Fiq)). 
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Figure 9: The function r = F{q) in the proof in section 3.5 representing the lower branch for (? < of the 
the relation (3.19[3). It is increasing as a function of /x and convex if /x = const. 7 and 7 is sufhciently small. 



One can verify the following for (q, r) on the graph of F: g2{T) > 0, gi{q) is bounded from below; 
g2{T) is negative and bounded from above; if q- < 0, then g'i{q) / g'2{'T) is positive and bounded 
from below (in fact this ratio reaches 00 at g = g_, where g'2{'T) = 0). As we have seen, we can 
guarantee (/_ < by making /i = const. 7, and then, if 7 is sufficiently small, we obtain F"{q) > 0. 

4 Stability of Harmonic Solutions 



To analyze the stability of the harmonic scattering solutions of the form (2.2), we first project 



the system to the resonator. This results in an equation for z alone exhibiting dissipation in the 
form of a delayed response coming from the coupling to the Schrodinger string with a point-mass 
defect. We then linearize about harmonic solutions Ze^**^*. Linear stability analysis is carried out 
by analyzing the determinant D{s) of this linear system in the Laplace-transform variable s in 
the regime of small 7 and fi. Any zero of D{s) in the right half plane indicates linear instability, 
whereas all zeroes being in the left half plane indicates linear stability. 

The responses |Z/ Jp are obtained from the roots of the polynomial f{P) ( 2.7|2.8 ), and it is the 
roots P themselves that appear in the expression for D{s). When 7 and fj, are small and u — Eq 
is bounded from below, one root of / is very close to zero, and when there are three roots, two 
of them are very close to 1 (Fig. [lOJ top). For the highest and lowest responses, the real parts of 
all zeroes of D{s) are shown to have a very small negative real part asymptotically as 7 — )• and 



/x ^ 0. A precise statement in terms of D{s) is made in section 4.4, and the result is the following 
theorem. 

Theorem 2. 

1. The harmonic response corresponding to the lowest root of the polynomial f{P) is linearly stable 
for uj — Eq > if ^ and fi are sufficiently small. 

2. When the system admits three distinct harmonic scattering solutions (necessarily ui — Eq > 0) 
and if 7 and /x are sufficiently small, then 

a. the solution corresponding to the middle root is linearly unstable. 

b. the solution corresponding to the highest root is linearly stable if oj is large enough and 
linearly unstable if Eq < uj < 1/2. 

Part (2b) has an intriguing consequence. If 7 and /x are small enough and chosen such that the 
ratio /x/7^ is large enough, there is a single frequency interval (wi, 00) of triple harmonic solutions 
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for which the highest response is unstable for small frequencies and stable for large frequencies. 
Thus a transition from instability to stability occurs at some frequency. 

Let us assume that, for i < 0, the system is in a harmonic scattering state {uh{x, t), yhit), Zh{t)) 



of the form (2.2). The state is perturbed for t > by forcing the resonator by a small-amplitude, 



temporally localized function e(t). Make the following substitutions in the main system (1.1): 



u[x, t) = Uh{x, t) + v{x, t), 
y{t) = yh{t)+7]{t), 
z{t) = Zh{t) + C{t), 

The deviation {v{x,t),r](t),C{t)) of the solution from the harmonic one vanishes at t = 0, as does 
e(t), and it satisfies the system 

ivt + Vxx = for X / 0, (4.36) 

iV = -fC-{vx{0^,t)-Vx{0-,t)) with r]{t) = v{0,t) , (4.37) 

iC = EoC + 77? + X{2\zhW + zU + 2zh\C\^ + ZhC^ + ICI'C) + e{t). (4.38) 



In addition, we impose on v an outgoing condition as |x| — )• (X), discussed below; see (4.42). This 
outgoing condition postulates decay of the v{x,t) as |x| — )• oo. It is symmetric in x because (i) 
V starts at rest {v{x,0) = 0), that is, u is a pure harmonic solution for t < and (ii) the value 
f (0, t) together with the requirement of decay determines v{x, t) for both x — )• oo and x — )• — oo. 

4.1 The outgoing condition 

The outgoing condition is understood through consideration of an auxiliary problem on the half-line 
X > without forcing and with a free endpoint at x = 0, 

ia = —Uxx, X > 0, i > 0, 

a(x, 0) = 0, X > 0, (initially at rest) (4.39) 

a(x, t) — ;■ as X — >• oo , t > 0. (decaying at oo) 

The Laplace transform of this system is isa = —hxx with a(s,x) — >• as x ^ oo. The solution 
satisfies hx = i^''^yj~sa, with arg(i^'^) = 37r/4, branch cut of ^ on the negative half line, and 
Re -y/s > 0, or 

a(x,s) = a(0,s)e*'^'v^^ (x > 0) (4.40) 

The inverse Laplace transform gives the general solution in terms of the Laplace-transformed value 
of a(0,t), 

Joo+O 

o(x,t) = — j a(0,s)e*'^'^^e^*ds. 
2m J 

— ioo+O 

The part of the integral along s = —iuj + (w > 0) is a superposition of radiating (outward 
traveling) waves, and the part along s = iw + (a; > 0) is a superposition of spatially evanescent 
fields. 

An analogous argument gives the outgoing condition for a function 6(x, t) defined for x < 0: 

6(x,s) = 6(0,s)e-*'^'v^% (x<0) (4.41) 

which is equivalently expressed as bx = —i^''^y/sb. 

20 



The value of v{0,t) connects v{x,t) to the left of the defect (x < 0) continuously with v{x,t) 



to the right of the defect {x > 0), that is, one puts a{0,s) = v{0,s) = 6(0, s) in (4.40) and (4.41). 



Thus the outgoing condition for the perturbation v is expressed in the Laplace variable by 

^ I —i^^'^y/s'v{x,s) for X < 0, /, , n 

Vx{x,s) = < I (outgoing condition) (4.42) 

I I' ysvyx^s) tor x > 0, 

in which arg(i'^/^) = 37r/4 and v{x, s) is continuous in x. In fact, v{x, s) is completely determined 
by t)(0,s). 

4.2 Reduction of the system to the resonator 



Because of (4.42 ), v is spatially symmetric, that is, v{x, t) = v{—x, t), and the jump in its derivative 
at X = can be expressed in a simple way: 



Vx{0+,t)-vx{0-,t) = 2vx{0+,t). 



Through the outgoing condition (4.42), this is expressed in the Laplace variable by 

C[v,iO+,t)-v,{0-,t)] =2^3/2, 



s rj . 



Equation (4.37) now yields a relation between rj and C, 

r?(t)=7£-i[5(s)C] (0=7(5*0(0, 
in which 



9{s) 



s + 2-\/is 



(4.43) 

(4.44) 

(4.45) 
(4.46) 



with the branch cut for J on the negative real half-axis and -y/r > for r > 0. The relation (4.45) 
allows one to project the system onto the resonator by considering equation (4.38) for a single 
function (^, 



K = EoC + i\g * C) + X{2\zh\^C + 4C + 2z^|CP + z^e + ICPC) + e{t) 



(4.47) 



The real part of the function g{s) is positive for s in the right half plane, which is a condition for 
power dissipation for a linear system (A = here) discussed in [5]. 

4.3 Linearization about a harmonic solution 



Equation (4.47) is linearized by eliminating the quadratic and cubic terms in (^ and replacing ({t) 
with the solution ^(i) of the resulting linear equation. It is convenient to remove the oscillatory 
factor e~*'^* and deal with the field ^(t) = ^(t)e*'^*. Keeping in mind that z^ = Ze"*"^*, one arrives 
at the following equation for ^■. 

iiP = {Eq - uj)^ + 7^^ * V' + A(2|Z| V + Z'^i)) + e{t)e'^^ , (4.48) 

in which p[t) = g{t)e^'^*. In the Laplace variable, this becomes 

{is + oj-Eo - T^p - 2A|Z|2)v; - XZ^i = e|,_i^, (4.49) 

{-is + u-Eo- j'^f) - 2A|Z|2)^ -XZ^ = #|,+i^. (4.50) 
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The second equation is obtained by conjugating the first, replacing s with s, and then using the 
rule f{s) = f{s) for Re (s) > 0. All quantities are analytic in s within their domains of definition. 
The determinant of this system is 

D{s) = 3X^\Z\'^ - 4A|Zp(a; - Eq - 7^ (Rep) ) + {-is + lo - Eq - -f'^p){is + lo - Eq - 7^^), (4.51) 

in which 

OJ + IS + ll\/ij: + IS 

P{s) = : 1 , ^ . (4.53) 

U — IS — l%\JbJ — IS 

The branch cut in the argument of the square root in the denominator of p is taken to be the 
negative imaginary axis and the branch is defined by \f\ = 1; this imparts a branch cut in the s 
variable along the half-line {s = si + iw, s\ < 0} in the left half plane. Enforcing the rule that 
p{s) = p{s) for Re (s) > dictates that the branch cut for the argument of the square root in 
p is the positive imaginary axis with ^/l = 1; this imparts a branch cut in s along the half- line 
{s = si—ibj, si < 0}. Thus D{si+is2) has two branch cuts along the half-lines {s = si^ioj, si < 0}. 
With these stipulations of the square roots, the denominator of p vanishes at the single point 
s = iuj and the denominator of p vanishes at the single point s = —ico. 

4.4 Stability analysis 

This section is dedicated to the proof of Theorem [2| which is stated in the proposition below in 
terms of the roots of D{s). 

The linear stability of the system about a scattering solution u^, that is, whether ^(t) grows 
or decays as t — )• cxd, depends on the roots of D(s). Any root in the right half s-plane indicates 
exponential growth, and all roots being in the left half plane indicates decay. In terms of the 
quantity P = fiTZ/a = X\Z\'^ /a, with a = uo — Eq — 'y'^/{uj + 4), D has the form 

D = 3P2_4P- (oj-Eo-^ip + p)] + ir{uj - Eq - is - -f^p){io - Eo + is - -f'^p). (4.54) 
a \ 2 J ^ 

Thus D{s^ depends explicitly on the parameters w, P, and 7^ (as well as Eq). The value of 
P is related to u^ through the correspondence between harmonic solutions and real roots of the 



polynomial /(P) := P(P — 1) +aP — /3 (2.8), which depends parametrically on 7, /x, and uj. When 
0. and /3 are small, the smallest root P\ is nearly zero and, in the case of three roots, the other two 
Pi and P2 are nearly 1. 

Proposition 3. 



1. In the expression (4.54) for D, let P he equal to the smallest root of f{P). If uj — Eq > is 
bounded from below and'j and fj, are sufficiently small, then all zeroes of D{s) have (small) negative 
real part. 

2. Suppose that f{P) has three roots (necessarily u > Eq). 



a. In (4.54), let P be set to the intermediate root. If j and /x are sufficiently small, then D{s) 



has a root with positive real part. 
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b. In (4.54), let P he set to the largest root. If the conditions 

1 

{u"^ - l){uj - E^f >u? 
47V 



(a;-^o)3(^ + 4) 
II OJ — Eq 



< Il-Pl <1 



72 Uj{uj + 4:) 

are satisfied and ifj and fx are sufficiently small, all zeroes of D{s) have (small) negative real part. 
If Eq < uj < 1/2, then D{s) has a zero with (small) positive real part. 

Let us simplify notation by putting 

p = u! + A, a = uj - Eq - -f'^/p , 



a 



474 



Loa^p"^ 



/? 



47V 
a^p 



P=l + Q. 



Case P < 1. 

We assume that uj — Eq > is bounded from below and let 7 and p tend to zero. In this regime, 
a r^ u! — Eq, so 1/a = 0{1). The definitions ( |2.9| ) of a and (3 show that these quantities vanish as 
7, // — )• 0, and thus the smallest root Pi of /(P) := P{P — 1)^ + aP — /3 is asymptotic to /3: 



P ^ 13 



a^p 



Oh'p). 



(4.55) 



Case P <^ 1 and p, p bounded. The first two terms of (4.54) vanish in this regime, and thus 
the third term also vanishes, yielding two cases. 



{cv- Eq 



is 



72^)^0 



or 



UJ 



Eo + is - 7^p ) 



0. 



Putting s = si + is2 in the first case yields 



S2 ~ — (cj — Eq) < . 



si 



-7 Imp. 



(4.56) 



(4.57) 
(4.58) 



The sign of —Imp is the same as that of — Re (si + 2-y/a7^— JsT^F^). This latter expression is 
asymptotic to —Re (2-^li; + S2) ~ —2y/Eo. This quantity is negative by the declaration of the 



square root in the definition of g (4.46). Thus si is asymptotically negative. A similar argument 
for the second of the cases (4.56) shows that S2 ^ oJ — Eq > and that si ~ 7^ Imp < 0. Thus the 
roots s oi D{s) = are in the left half plane. 



Case P <^ 1 and p or p unbounded. As we have mentioned at the end of section 4.3, the 
denominator of p vanishes only at iuj and that of p only at —iuj. Thus, if one of these quantities 
is unbounded, the other remains bounded. It suffices to analyze the case of p being unbounded, as 

p{.s) = p{s) and D{s) = D{s). 
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If Re{s) > 0, then \/uj + is is in the upper half plane. The square root in the definition (4.52) 
of p takes arguments in the upper half plane into the first quadrant, and thus 



(4.59) 



Re{s) > =^ Imyjuj + is > . 

We will prove that this condition is asymptotically inconsistent with D{s) = 0. 

The assumption that p — )• oo implies s — )• iu. Applying this to the three terms of D{s) in (4.54) 
gives 



4, ,2 



D 



487> 



:i+o(i)) 



167^/x 



7 



.V v'-y''i<'+°™+ 



+ 



2uj- Ei 



a^ 



(i^o(l + o(l))-7'p)(l + o(l)) 



7^P 



O" 



^(2u;-£;o)(l + 0(1)) + ^(1 + 0(1)) 



a-p 



+ ^(2a;- ^0)^0(1 + 0(1)) + 



(T^ 



I67V.-, ,.^^ 487V^, 

^^(l + o(l)) + ^^(l + o(l)) 
ap a^p^ 



l^P 



l-{2u;-Eo){l + o{l))] + :^(2a;-i?o)(l + o(l)). 



a^ 



Setting D = provides the asymptotic relation 



En 



P ^ ^2 ■ 

T 



On the other hand, the definition (4.52) of p with u; + is — >• provides the relation 

1 

P 



2i\/w + is 
Combining these two asymptotic expressions for p yields 

„2 



\^^ 



IS ~ —I 



7 
2En 



which, in view of (4.59), is inconsistent with Re (s) > 0. 
Case P ~ 1. 

To see the asymptotics of the other roots, write f{P) = in terms of Q = P — 1: 

47' /^ 



f{P) = 
This implies the asymptotic 



^ P \7 



ujpj 



Q'-'4(^--) -hen ^«M«1. 
o' P V7 ^PJ 



(j'^p 



(4.60) 



Thus there are two values of order 7^ that Q can take on, one negative and one positive, under two 
conditions: 

^-^ < IQI < 1 and ^> — 



a-'p 



7^ ujp 
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Q~±C772 (C>0). 



(4.61) 



Equality in place of ">" in the condition /i/7^ > a /{up) is achieved asymptotically at the 0J1-0J2 
bifurcation. In general, for a fixed asymptotic ratio /x/7^, the inequality of satisfied if uj is large 
enough or a is small enough. 



The expression (4.54) for D{s) can be written as 



D = 2Q-^ + ^{p + p) 
ap a 

-^Q-(--lip + p))+^ 

a \p 2 J a'^p'^ 



^H [P-P) + 



4 4 

7 - 7 - 

-^{p + p) + -^pp- 



(4.62) 



Case P r-^ 1 and p, p bounded. These assumptions imply 



ap a 



4 + ^(^-p) + 0(7^). 



(4.63) 



Let us assume Q ~ ±6*7^ with C > from conditions (4.61). Setting D to zero, one obtains 
\s\ <C 1, and expanding p and p in s gives 



p + p 
p — p 



2s 



p pco^/'^ 



puj 



1/2 



+ 0{\s\ 
0{\s\). 



With these expressions, the equation D = becomes 
-2^=0 (1 + 0(7")): 



o + 2,^^^<^) +Oh*) 



pu 



3/2 



Because Q ~ ±6*7^ with C > 0, the 0(7^) on the right-hand side may be absorbed into the 0(7^ 
on the left-hand side, 

s = V-20-2Q -27^ 



a[u) 



puj 



3/2 



+ 0(7^), V-2<t2Q ~c7 (c/0). 



(4.64) 



The negative root Q ~ —C^"^ corresponds to the middle root P2 of f{P), or the intermediate 
response of the resonator depicted by the middle branch of the amplitude vs. frequency graph at 
the top of Fig. I5J Thus D{s) has a zero in the right half plane for sufficiently small 7. 

The positive root Q r^ Cj'^ corresponds to the largest root P3 of f{P), or the highest response 
of the resonator depicted by the top branch of the amplitude vs. frequency graph. The first term 



of (4.64) is imaginary, so the second term, which is real, determines the sign of the real part of s. 



If Eq < 1/2, then one has Re (s) > for 00 < 1/2; otherwise. Re (s) < 0. 

Case P ~ 1 and p or p unbounded. Again, it suffices to analyze the case of p being unbounded. 



Let us suppose that p is unbounded and p is bounded. The balance of dominant terms in (4.62) 
yields 



7^ 
a 



(l + is) 



G^ 



(4.65) 



Formula (4.52) with p — )• 00 gives 



VjJ , p 



2i\Juj + is 
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Using the second of these in (4.65) gives 



VoTT 



IS ~ 



7^cr(l + is) 



2is^ 



and then using s ^ iuj in the right-hand side yields 



VoTT 



IS 



ij'^a^ca - 1) 



2^2 



as long as (jj ^ 1. At a; = 1, (1 + is) = o(l). In any case, we obtain uj + is = 0(7^), so that 

s2 = _a;2 + 0(74). 



Each of the terms is and s^ appears only once explicitly in ( |4.62 ), and they may be replaced by —uj 
and — a;2 committing an error of only 0(7^). Let us also introduce the proper scaling of Q from 



the condition (4.61), namely Q = 7 Q 



v2n. 



ap a a-^ a 



(4.66) 



-4Q f^ - ^{p + p)] + 4^ - ^{P + P) + \pp + 0{^^) . (4.67) 



Passing all terms of order 7^ into the error (recall that p = 0(1)) and rearranging terms to isolate 
the quantity of interest 7^^ yields 



D 



1 " 
IP 



1 - w + 7^ 2Q 



ap 



UJ 



a^ 



■J , ~ 2 1 , w „ 

7^ 2Q + -p + -p 

ap a a 



0(7^) . (4.68) 



Now setting D = gives 



l^P 



a(w-l) 



W2 



2Q 



2 1. w., ,^, 
— + -P+-P +0(7 
op a a 



l + ^f2Q + ^?^^)+0(/^ 



a;-l 



o-p 



(4.69) 



The real and imaginary parts of this quantity are 



-A := Rej'^p 



J B := Im7 p 



UJ 



+ 0(7'), 



a{uj-l) 
j'^Lo'^lm.p / {u> + l)a 1 



+ 0(7' 



(t{uj — 1) \ a;2 (w — l)cr^ 

in which jj is evaluated asymptotically using u + is = 7^^, where ^ is bounded: 

1 1 

p 



V2w - 7^^ {-2i + v^2w - 74^ ) \/2^(-2i + ^/2^) 
The quantity 72^ can be evaluated alternatively directly from its definition (|4.52): 



1 



V^7TI^(2i + \/Z7T^) 272^1 (1 - i7^\/?) 

272 VfV 2 4 + 



(4.70) 
(4.71) 

(4.72) 



-i ^ 1 ^ ^7^\/? ^ 



272^1 4 8 
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which results in the asymptotic 



7^p= -^ + ^ + 0(/). (4.73) 



By equating expressions (4.69) and (4.73) for 7 p, we obtain 

in which A and i? are real and differ from A and B by order 0(7^). If uj is sufficiently large, 
A and -B are positive and bounded from below in magnitude. Specifically, it is sufficient that 
(w^ — l){uj — Eq)"^ > io'^ . This yields Imp > asymptotically. 

Finally, 

IS = —UJ + 74 = ~^ + 



4 (74^2 + ^^2)2' 
which results in Im s ~ w and 

This result places the zeroes of D{s) asymptotically in the left half plane, regardless of the sign 
ofQ. 

5 Discussion of continuum-oscillator models 

Simple continuum-oscillator systems serve a vital role in elucidating fundamental principles and 
phenomena in physics. Horace Lamb, interested in how disturbances in a body subside due to the 
transmission of energy into an infinite ambient medium devised what is now known as the Lamb 
model as the simplest expression of this phenomenon |8]. He showed that, if a harmonic oscillator 
is attached to an infinite string whose displacement is governed by the wave equation, the oscillator 
obeys the equation of the usual instantaneously damped harmonic oscillator. In other words, the 
energy loss in an oscillator due to instantaneous friction can be perfectly conceived as the radiation 
of energy into an infinite string. The coupled system of the oscillator and the string together is a 
conservative extension of the lossy system consisting of the damped oscillator alone. It is in fact 
the minimal conservative extension of the damped oscillator, and it is unique up to isomorphism, 
as shown by Figotin and Schenker [5]. The string acts as a system of "hidden variables" from the 
point of view of an observer who is able to make measurements only of the motion of the oscillator. 



The linear version (A = 0) of our model (1.1) is designed specifically to allow the oscillator to 



be completely detached from the string in the zero-coupling limit (7 = 0) without severing the 
string. The case 7 = corresponds to a resonator decoupled from a system of "hidden" variables 
that itself exhibits non-resonant scattering in a line by a point-mass defect resulting in a simple 
but nontrivial transmission coefficient. A very small coupling parameter 7^1 corresponds to a 
small perturbation of a decoupled system with an embedded eigenvalue Eq and results in a sharp 
resonance. 
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Figure 10: Top. Each of the three roots P^, i = 1,2,3 of the polynomial f{P) = P{P ~ l)"^ + aP - 13 
(2.8), all of which are positive, corresponds to a harmonic solution of the system ( |1.1[ ). The response of 
resonator is related to P by equation (2.7). With the parameters 7 = 0.05, ^ = 47^/^23 ~ 0.000429, and 
Eq = 4, and w = 1.5 0^23 ~ 14.5, Pi is very close to 0, whereas ^2,3 are very close to 1. (There is only one 
transition frequency uii for these parameters, i.e., fj, is above the point B23 in Fig. l4J) Bottom. For each 
root Pi , the zero sets of the real and imaginary parts of D{s) (4.54) in the complex s-plane are shown from 



left to right in increasing order of Pi. These graphs only show the zeros of D{s) that are away from ±iu}. 
For 7 sufficiently small, the imaginary part of D vanishes on the real axis and on an almost vertical curve 
approximately connecting the points —iuj and zw, whereas the real part vanishes on a hyperbola-like curve. 
The real part of this curve is small and negative, in agreement with the asymptotic calculation in section [44] 
when p OT p is bounded. 
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Because our transmission line is governed by a Schrodinger equation, it exhibits dispersion, 
which is experienced by an observer in the oscillator as a delayed response, or a non-instantaneous 
friction, assuming that there is no forcing originating at points along the length of the string: 

/■oo 

iCit) = Eom+7^ g{t')C{t-t')dt' + e{t), (A = 0) 
Jo 



9{s) 



-I 



S + 2^/is 



The function g{t) obeys a power-dissipation condition, described in [5j, which in the Laplace variable 
is expressed by the condition that g{s) has a positive real part when Re (s) > 0. 

The function e(t) in the above equation is a spatially and temporally localized perturbation 
of a harmonic oscillation, that produces a deviation C,{t) in the state of the oscillator. When 
the frequency of oscillation vanishes, the problem becomes that of the dissipation of finite-energy 
disturbances of a system initially at equilibrium. In this case, nonlinearities of a general form 
have been analyzed by Komech [7j when the string's motion is governed by the wave equation 
and the resonator is attached as in the Lamb model. The system exhibits transitions between 
stationary energies of the nonlinear potential in the resonator, which resemble transitions between 
energy states in atoms. In Komech's model, the energy is related to the height of the string rather 
than a frequency of oscillation. A positive cubic nonlinearity has only one stationary point and all 
disturbances decay to zero. 

At nonzero frequencies, cubic nonlinearity becomes interesting and the focus of study turns to 
the steady-state behavior of a nonlinear scatterer subject to a monochromatic harmonic forcing 
originating from a source far away. In particular, one wants to understand how these steady oscil- 
latory motions respond to finite-energy perturbations. Our choice of a Schrodinger equation for the 
string was based on the form of the nonlinearity A|zpz that is natural for this equation, and which 
admits periodic solutions that are purely harmonic and mathematically tractable. Perturbation 
about the harmonic motion Ze"**^* of the oscillator results in the equation 



oo 

iu)t 



ill, = {Eo - uj)^) + 7^ / p{t') il}{t - t') dt' + A(2|Z| V + ^V) + e{t)e 
Jo 



(4.48) for the linearized perturbed oscillation ^(t)e~*'^* of the nonlinear system. Energy loss comes 
from the delayed response term p{s) = g{s — iuj), which also satisfies the dissipation condition that 
Re (s) > =^ Re (p(s)) > 0. The solutions Z depend in a complex way on 7 and A, and stability 
analysis when these parameters are small is delicate, as demonstrated in section |4j 

The idea of projecting a conservative oscillatory system onto a lossy subsystem is an insightful 
one and has been discussed from different points of view in the works mentioned above. Figotin 
and Schenker [5] view the whole string-oscillator system as a conservative extension of a dissipative 
subsystem of "observable variables" (the oscillator). The string realizes in a structurally unique 
way a space of "hidden variables" that are responsible for the loss of energy measured by an 
observer confined to the oscillator. The form of the energy dissipation (as the function g(t) above 
in the linear case) observed in the oscillator is sufficient to determine the space of hidden motions 
responsible for the dissipation. These ideas are brought to bear on the lossy Maxwell system in 
electromagnetics by these authors and others, such as Tip [18]. Komech describes the projected 
system as an irreversible description of a larger reversible one. This point of view had been advanced 
previously by Keller and Bonilla [6] as an illustration of how irreversible processes may be derived 
from reversible ones and was motivated by the question of whether macroscopic physical processes 
can be deduced from classical mechanics. 
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To an observer of our nonlinear system from the site of the resonator, the couphng to the string 
is felt as a combination of input energy and energy loss to damping. These two energies balance 
out (in the sense of time averages), when the system is in a harmonic steady state. The balance is 
disturbed when the system is perturbed from steady state. To express this perspective in precise 
terms, we project the system onto the resonator. Write the equation for tp above in terms of the 
actual perturbation (^(t) = ip(t)e~^^* of a harmonic solution Zh = Ze~"^* and reinstate nonlinear 



terms from (4.47): 



/•CO 

K = EoC + 7^ / 9{t') Cit - t') dt' + A {2\Z\^C + Z^e-^^'^'C) + 
Jo 

+ A (2Ze~^'^*|Cp + Ze^'^'C'^ + |CpC) + e{t). 

Both the external forcing e and the steady-state field Zh = Ze'^'^^, induced by the incident field 
j^i{kx-ujt) ^ affect the dynamics of the system. The field z^ = Ze~^^^ depends on 7, //, and uj 



through the roots of the cubic polynomial /(P) (2.8). If e{t) is taken to be eo'5(i), where 5(t) is 



a unit impulse at t = 0, we can consider this equation for t > with e = and a nonzero initial 
condition. 

Importantly, introducing external forcing and damping through coupling to the string leads to 
harmonic solutions that can still be calculated explicitly. This is not possible if external forcing 
and instantaneous damping are introduced directly as in the much studied Dufhng oscillator 

z + az^ Eqz + 6/ = Fo cos(a;t) . (5.74) 

In spite of the similarities of cubic nonlinearity, harmonic forcing and intervals of triple solutions, 
these solutions are only approximately periodic in the Duffing case, when the parameters a, h and 
Fq are small. 
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